Exact temporal evolution of the two-species Bose-Einstein condensates 



o 

(N 



oo 



c/2 

^— > ' 

c : 

ctf . 
=5 
CT 
-4-* ! 



i 

C 

O 
o 



> 

(N 



m 
O 
m 



X 



Cong Zhang, Zhi-Hai Zhang, and Shi-Jie Yang* 
Department of Physics, Beijing Normal University, Beijing 100875, China 

We construct exact stationary solutions to the one-dimensional coupled Gross-Pitaevskii equations 
for the two-species Bose-Einstein condensates with equal intraspecies and interspecies interaction 
constants. Three types of complex solutions as well as their soliton limits are derived. By making 
use of the SU(2) unitary symmetry, we further obtain analytical time-evolving solutions. These 
solutions exhibit spatiotemporal periodicity. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) in weakly inter- 
acting atomic gases have offered a practical means of 
studying nonlinear behavior by using the matter waves. 
Among the macroscopic excitations, solitons and vor- 
tices are the most common and are extensively stud- 
ied both theoretically and experimentally. The repulsive 
BEC has received most experimental attentions while 
the attractive BEC is believed to only be stable in one- 
dimensional (ID) systems p], H|. The Gross-Pitaevskii 
equation (CPE), which describes the mean field dynam- 
ics of a dilute BEC, can be reduced to a quasi-lD nonlin- 
ear Schrddinger equation (NLSE), with a limit that the 
transverse dimensions of the condensate are on the order 
of its healing length and the longitudinal dimension is 
much longer than its transverse ones. 

The observation of solitons in BECs brings great pro- 
gresses in the study of the exactly integrable nonlin- 
ear systems which have a variety of applications in 
fiber optics as well as other fields Dark solitons 

are expected only to exist in repulsive interactions (s- 
wave scattering length a > 0), whereas bright soli- 
tons exist in attractive interactions (a < 0). The 
ID NLSE is exactly solvable by the inverse scattering 
transform 0, and has a countably infinite number of 
conserved quantities 0, [l(J. For the single-component 
BEC, the stationary solutions have been deduced analyt- 
ically under box or periodic boundary conditions (a. 1 1 lfl . 

A very relevant generalization of this class of physical 
systems is the coupled multiple-species BECs. Various 
soliton complexes including bound dark-dark jl2j. dark- 
bright 13J, dark-antidark, dark-gray, bright- antidark, 
and bright-gray [3| were predicted. In Ref-|15|. the au- 
thors constructed solutions of the coupled NLSEs with 
additional linear Rabi coupling to the problem where this 
linear coupling is absent. In this paper we seek exact 
time-evolving solutions to the ID coupled GPEs which 
describe the two-species BECs. For equal intraspecies 
and interspecies coupling constants, the equations of mo- 
tion are an integrable Manakov system which is applica- 



ble to the BECs consisting of different hyperfine states of 
87 Rb. In this case, the Hamiltonian has the SU (2) sym- 
metry. The normalized GPEs for the mean-field order 
parameter ^ = (ipi, ^2) T in the uniform external poten- 
tial read, 
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fi£*(M)+7l*(*>t)l a *(M), (1) 



where 7 is the intraspecies and interspecies coupling con- 
stant. 

The stationary equations are obtained by substituting 

tpi{x,t) = ipi(x) exp(— ifXit) (i — 1,2) as 

IHhix) = -\d 2 Jj t {x)+ 1 { J2 \^Ax)\ 2 )Mx). (2) 

.3=1,2 

It should be noted that the solutions to Eqs.Q are not 
really "stationary" because the time-dependent phase 
factor in each component is different as /ii ^ fX2- So 
there exists a Lamor procession around the z-axis in the 
pseudo-spin space. However, the density distribution of 
each component is time- invariant. Hence we briefly call 
them the stationary states. The periodic boundary con- 
ditions 

^(l)=Vi(0), $(1)=$(0), (3) 
are applied and the wavefunctions are normalized to 

(\Mx)\ 2 + \Mx)\ 2 )dx = l. (4) 

We first present three types of complex-real solutions 
to the stationary equations ([2]) and examine the corre- 
sponding soliton limits. By applying a unitary trans- 
formation to Eq.|T]), we can obtain exact time-evolving 
solutions. 



II. TYPE A 

We first consider the following form of stationary solu- 
tions to Eq.@ by the combination of function sn(fcx, to) 
and cn(/cx, to) as 
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^(x) = Dsn(kx,m), 



(5) 
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where fix) — y/A + Ben 2 (kx,m) with A, B, D the real 
constants, sn and cn are the Jacobian elliptical func- 
tions with modulus < to < 1. The period [HI, [13] 
k = 4jK(m) with Kim) is the complete elliptic inte- 
gral of the first kind. In the following we always choose 
the number of periods j = 2. By substituting the 
solution ([5J into the stationary Eqs. ([2]) and making 
use of the identity between the Jacobian elliptical func- 
tions sn 2 + cn 2 = 1, we obtain the following decoupled 
equations [H, FToj ] 



MiV'iO) = --kdltpiix) +7i|?/>i(x)| 2 V>i(2;) 
fail)2{x) = -\dl^iix) + 72 1 ^2(2;) | 2 ^(x), 



(6) 



where /xi^ are the effective chemical potentials and 
71,2 the effective intraspecies coupling constants in each 
species, respectively. One has 



Mi =m-~{{A + B)^ 
fo = tX2--y{A + B), 



(7) 



The periodic boundary conditions ([3]) require 



(1) -6(0) = 2jir x n, 



(12) 



where n is an integer. The periodic condition for the 
phase 6(x) can be fulfilled by properly adjusting the 
modulus to of the Jacobian elliptic functions. Figure 
1 (c) and (d) show the distribution of the phase for 
to = 0.1 and to = 0.999, respectively. The Fig. 1(d) and 



and 




7i =7(1" f) 
72=7(1--Jr)- 



(8) 



The decoupled Eqs.© can be self-consistently solved 
to yield 



Mi = ±fc 2 (l - 2?n 2 ) - l^ 1 A + 1 (A + B)%- 



M2 



ifc 2 (l + to 2 ) +7(^4 + 5), 



(9) 



FIG. 1: Left two columns: temporal evolution of the density 
profiles for the state (|14[) . For the non-soliton cases (upper 
row), the parameters are m = 0.1, 7 = 100, /ii = 41170.48, 
fi 2 = 41171.53. For the soliton cases, m = 0.999, 7 = 0.01, 
[ii — 1322.01, and {12 = 1293.58. Right column: the phase 
profiles of the stationary state © 



and 



D = 



B. 
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(10) 



There are two cases according to the value of B > 
and B < 0. For B > 0, we have the effective couplings 
71 < and 72 > 0, independent of the sign of the real 
nonlinear interaction 7. For B < 0, the nonlinear inter- 
action 7 must be repulsive and the effective couplings of 
both species also should be repulsive (71 > 0,72 > 0). 

The phase 9(x) is determined by the imaginary part of 
the decoupled Eqs.© which is given by 



(x) = 



/ 2 (0 



(11) 



where a = ±(2/ii^l 2 - 2jiA 3 + fc 2 (l - m 2 )AB)i is an 
integral constant. 



(f) are the phase profiles 9(x) of stationary state "J. The 
monotonous gradient in the phase profile indicates that 
a supercurrent is carried by %p\. 

We can further obtain the exact time-evolving solu- 
tions to Heisenberg equations of motion (IT]) by making 
use of the SU(2) symmetry of the Hamiltonian. We con- 
sider an arbitrary pseudo-spin rotational transformation 
y'(x,t) = U(<f>)y(x,t) withlH 



U(<f>) 



cos <p sm < 
— sin <b cos < 



(13) 



where the time factor e _IMi * is attached to each compo- 
nent. Consequently, ty'{x,t) is also a solution to Eq.(fT]) 
which is explicitly written as 



ip' x ix,t) = exp(— ifj,2t)[yfA + Bcn 2 ikx,m)e ie ^ cos</>exp(— ifxt) + 
ip' 2 (x,t) = exp(— ifi2t)[— \J A + Bcn 2 ikx, m)e l9 ^ sin^>exp(— i/it) 4 

I 



Dsn(kx, to) sin <p 
Dsnikx, to) cos ( 



(14) 



where fj, = — fj,2 is the temporal evolving frequency. Obviously, the density distribution possesses the spa- 
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tiotemporal periodicity. 

Figure 1 shows the temporal density evolution of the 
state *&'(x,t) (fT4|) . Hereafter we always choose the ro- 
tational angle 4> = n/8. Fig.l(a-c) are for m = 0.1 and 
the corresponding parameters are 7 = 100, A = 101.38, 
B = 309.5, D 2 = 309.6, and n = 1. Fig.l(d-f) are 
for m ~ 0.999 which correspond to the soliton limit. 
The relevant parameters are 7 — 0.01, A = 3.3 x 10 3 , 
B = -3.2 x 10 3 , D 2 = 1.26 x 10 5 , and n = 3. It displays 
a dark-gray solitonic complex. 



The relevant parameters are 



Mi = ifc 2 (l - 2m 2 ) - §71,4 - 1 A^ f 



and 



11-2 



= \k 2 {l-2m 2 ) + 1 A, 



(18) 



III. TYPE B 



D = — 



-B. 



(19) 



By the same way, we construct the following form of 
complex solutions 



^{x) = Dcn(kx,m), 



(15) 



where f(x) = \J A + _Bcn 2 (kx,m) with A : B, D the real 
constants. The effective chemical potentials jli^ and ef- 
fective coupling constant 7 12 are obtained as 



and 



jSi = hi + -fA^- 

P-2 = - 7 A 



71=7(1 + ^) 

72 = 7(1 + Jr)- 



(16) 



(17) 



The phase 9(x) is also expressed by Eq.(fTTj). Figure 2(c) 
and (f) show the phase profiles of the state (TTSl) for m — 
0.1 and m = 0.999, respectively. It indicates a super 
current in the condensate. 

There are also two cases for B > and B < 0. For 
B > 0, we have 7 < and the effective couplings of the 
two species should be attractive, 71 < and 72 < 0. For 
B < 0, we note that the effective coupling in the first 
species is repulsive (71 > 0) while the effective coupling 
in the second species is attractive (72 < 0), independent 
of the sign of the real nonlinear interaction 7. 

After pseudo-spin rotational transformation, we again 
work out the time-evolving solution to the Heisenberg 
equations of motion (TTJ for the two species as 



ipi (x, t) = exp(— i/i 2 i) [yC4 + Bcn 2 (kx, m)e 10 ^ cos tj> exp(-ifit) + Dc\i(kx, m) sin cj) 
^(x, t) — exp(— i(i2t)[— \J A + Bcn 2 (kx, ni)e l6lyX ^ sin</>exp(— i/it) + Dcn(kx, m) cos< 



(20) 



r 



where fi = /ii — [12- 

Figure 2 shows the temporal density evolution of the 
state ^' (|2H)) . Fig. 2(a-c) correspond to non-soliton 
cases m = 0.1 with the attractive coupling 7 = 0.1, 
A = 517.5, B = -287.26, D 2 = 271.3911 and n = 3, 
while Fig. 2(d-f) are for the soliton limits m — 0.999 with 
repulsive coupling 7 = —0.1, A = 1031.5, B = —865, 
D 2 = 1.3774 x 10 4 and n — 4. Under this circumstances, 
we obtain a bright-bright solitonic complex in the con- 
densate. 



IV. TYPE C 

We can construct the third type of complex solutions 

(21) 



^(x) = f{x)e w ^ 
^{x) = Ddn(kx,m). 



The effective chemical potentials and effective coupling 
constants are respectively 



Mi = ^i-iD 2 {l + m 2 )- 1 M 
fi 2 =H2-lA- 1 B^ 1 



(22) 



and 



(23) 



71=7(1 + ^) 
72 =7(1 + ^7). 

The decoupled equations © are self-consistently 
solved to yield 

/ tlx = |fe 2 (l - 2m 2 ) - \j,A + 7^ 2 (1 + m 2 ) + ^A^^ 



fi2 = lk 2 (m 2 -2) + 1 A + 1 B^, 



and 



B + m D 



2n2 



m k 



21,2 



(25) 



XXX 




FIG. 2: Left two columns: temporal evolution of the density 
profiles for another state (|20|l . Here are the given parameters: 
(a,b,c) the non-soliton state shown in upper row corresponds 
to m = 0.1, 7 = 0.1, Hi = 130.90, (jl 2 = 129.51; (d,e,f) the 
soliton cases corresponds to m = 0.999, 7 = —0.1, fii = 
22.3771 and {12 = —747.29. Right column: the phase profiles 
of the stationary state (|15[) . 



The signs of the effective couplings in the two species 
are similar to those in the TYPE B solution. 



with the temporal evolving frequency fi — fix — M2- 

In Fig. 3, a graphical temporal evolution of the density 
profiles of this solution is displayed for 7 = —0.1. Fig.3(a- 
c) demonstrate case of m — 0.1 with the corresponding 
parameters are calculated as A — 555.5, B = —308.36, 
D 2 = 2.925 x 10 4 and n = 3. Figure 3(d-f) display the 
soliton limit with m — 0.999. We obtain the parameters 
A = 60.268, B = -57.92, D 2 = 1.2992 x 10 4 and n = 3, 
respectively. The condensate also forms a bright-bright 
solitonic complex. The phase profiles 0(x) of stationary 
state \1/ is displayed in Fig. 3(d) and (f). 

V. SUMMARY 

In summary, we obtained three types of stationary so- 
lution to the pseudo- spin- 1/2 BEC which possesses the 
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Through the pseudo-spin rotational transformation, 
we obtain the time-evolving solutions to the Eq.(p} as 



(26) 

I 

SU(2) symmetry. Exact time-evolving solutions are con- 
structed through a pseudo-spin rotational transforma- 
tion. The states exhibit periodicity in both space and 
time. In the limit of m — > 1, we obtained the solitonic 
complexes for the two species condensate. 
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